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Stacks of intrinsic Josephson junctions in Bi2Sr2CaCu208+5 emit intense and coherent terahertz 
waves determined by the internal electromagnetic cavity resonance. We identify the excited trans¬ 
verse magnetic mode by observing the broadly tunable emissions from an identical nearly square 
stack and simulating the scattering spectrum. We employ a wedge-type interferometer to measure 
emitted integral power independently of the far-held pattern. The simulation results are in good 
agreement with observed resonance behaviors as a function of frequency. 

PACS numbers: 74.50.+r, 74.72.-h, 85.25.Cp 


I. INTRODUCTION 

Electromagnetic (EM) waves in the 0.3-10 THz fre¬ 
quency range are thought to have great potential in re¬ 
search and industr}^! thus, compact, solid-state, and 
continuous wave (CW) terahertz sources have been 
developed in the held of semiconductors and laser^. 
Presently, quantum cascade lasers, which emit radiation 
at frequencies greater than 1.45 THz and must be cooled 
to 37 kP, are of particular interest. Room-temperature 
terahertz sources based on resonant tunneling diodes can 
generate frequencies up to 1.31 THJ^but emission power 
greater than 1 mW is presently not possible. 

Intense, continuous, and broadly-tunable terahertz 
generation from high transition temperature (high-Tc) 
superconductor Bi2Sr2CaCu208+(5 (Bi-2212) has at¬ 
tracted the attention of experimenters and theoreti¬ 
cian^. DC voltage Vb applied to the stack of intrin¬ 
sic Josephson junctions (IJJs)P in single crystal Bi-2212 
leads to AC current and concomitant terahertz emission 
at the Josephson frequenc 3 pl 


where e, h, and N are the elementary charge, Planck’s 
constant, and the number of active IJJs, respectively. 

To generate intense terahertz waves, several tens of 
micron-sized mesa structures are milled out of Bi-2212. 
According to the electromagnetism for the microstrip an- 
tennsP, the transverse magnetic (TM) mode is excited at 
the resonance frequency: in the case of a rectangular 
mesa. 



where cq, e, w, and i are the speed of light in vacuum, 
the dielectric constant, and the width and length of the 


mesa, respectively. The emitted far-held waves are polar¬ 
ized with their E-held parallel to the plane that contains 
the wave vector of the standing EM wave in the mesa 
cavity. The two integers (m,p) in Eq. (2) correspond to 
the number of line nodes of the TM (m,p) mode stand¬ 
ing wave. In most cases using long rectangular mesas, 
the TM (1,0) mode with formation of the standing wave 
along the width of the mesa at /f q = co/2y^u; has been 
observecP^. The similar collective resonance-like be¬ 
havior along the length was found in laser microscop}^!!, 
which, however, is not direct evidence for the cavity res¬ 
onance. To our knowledge, there still remains a lack of 
consensus on the issue on the cavity resonance condition. 
Moreover, broadly tunable e missio ns obscure the role of 
the internal cavity resonanc d^^ l ^^ i To correctly identify 
the excited modes, the emitted integral wave has to be 
collected using lenses or mirrors with small f number be¬ 
cause the far-held emission patterns are affected by the 
excited modd^. 

A complete study of the cavity resonance will help to 
better understand the mechanisms of coherent emission 
from the IJJ stack, which will make possible to control 
the polarization of the emitted wave and may facilitate 
the design of CW terahertz sources for coherent commu¬ 
nication applications. In this study, we show that various 
cavity modes with different {m,p) can be excited using 
an identical stack of IJJs by changing the bias condition. 
We simulate the experimental data using the method of 
moments to calculate the resonance frequency and visu¬ 
alize the EM standing wave in the emitting mesa. 

II. EXPERIMENTAL SETUP 

We show experimental results from a nearly square 
mesa sample with w = 100 /im, i = 140 /im, and 
t = 1.3 um. The sample preparation method is described 
in Ref.“ All sample dimensions are measured using an 
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FIG. 1. Top view of the collection optical system. 


atomic force microscope. The horizontal dimensions (re, 
t) are averages at the top and bottom lengths. Figure 1 
shows the top view of experimental setup. The sample 
was mounted on a He-flow cryostat. The bath temper¬ 
ature Ti) was monitored using a thermometer attached 
to the sample holder. For measuring the current-volt age 
characteristics (IVC), we used a conventional electrical 
circuit for measuring the resistance, where a load re¬ 
sistance of 100 11, a standard resistor of 10 H, and the 
emitting sample were connected in series to the constant 
voltage source. During the I-V measurements, we also 
monitored the far-infrared emission power using a sil¬ 
icon composite bolometer. The bolometric signal was 
calibrated separately using a sub-millimeter power me¬ 
ter. For the lock-in detection, the emitted EM waves 
were optically chopped at 80 Hz. To collect the emitted 
integral waves, we used hemispherical silicon and plas¬ 
tic lenses and an off-axis parabolic mirror, as shown in 
Fig. 1. The collimation effect owing to the silicon lens 
is described in detail in the Appendix A. Note that the 
estimated solid angle of the detection window of 0.6 sr 
is more than 20 times larger than that in previous ex¬ 
periments; therefore, we can directly measure the total 
emission power from the sample. At the focal point of 
the parabolic mirror, we set a wedge-type interferometer 
system to measure the emission frequency The de- 
tails of the frequency measurement are described in the 
Appendix B. 


III. RESULTS AND DISCUSSION 


Figure 2(a) shows the I VC curve at = 30 K obtained 
using a cyclic bias scan. We subtracted the contact re¬ 
sistance owing to the two terminal measurements from 
the I-V data. In Fig. 2(b), the emission power on the 
outermost I VC curve, which was simultaneously moni¬ 
tored using the bolometer, is shown as a function of the 
current. The inset in Fig. 2(b) shows a scanning ion mi¬ 
croscopy image of the emitting sample. Four emission 
peaks are observed at 22.1, 15.0, 10.2, and 6.54 mA, as 
indicated by the arrows. Each emission peak corresponds 
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FIG. 2. (a) Gurrent-voltage characteristics at Tb = 30 K, 

obtained by the cyclic bias scan, (b) Emission power of the 
outermost IVG curve measured using the bolometer as a func¬ 
tion of I. The inset shows a scanning ion microscopy image 
of the emitting sample. 


to the excitation of the different cavity mode and is dis¬ 
cussed elsewhere. We scanned the bias current upward 
and downward several times to confirm the emission re¬ 
producibility and found that the emission peaks always 
appeared at the same bias points. 

The emission frequencies /(V) at = 30 K are shown 
in Eig. 3 together with the emission power. The er¬ 
ror bars for / are standard deviations from the sinu¬ 
soidal fitting in the interferometer analysis. The emis¬ 
sion power in the right axis is the same data as in 
Eig. 2(b). Three bias points at 0.92, 1.11, and 1.26 V 
correspond to the emission peaks at TM (0,2), (1,2), 
and (2,0), respectively. Note that / for the rightmost 
peak of TM (2,1) cannot be obtained properly because 
the bias point around 1.35 V is less stable to the bias 
scan. The dashed line represents the Josephson relation 
/ = fj = {2e/h)V/N, where V = 851 is obtained from 
the least squares fit to the experimental data. This N- 
value is consistent with the value of 850 that is separately 
obtained from the height of the stack. 

According to the Josephson relation, the bias voltage 
can be directly converted into the emission frequency. 
Eor instance, the application of 1.00 V voltage to the N = 
851 IJJ stack corresponds to a frequency of 0.568 THz. 
Hence, the plot of the emission power versus V can be as¬ 
sumed to be equivalent to the resonance spectrum. This 
assumption is supported by the fact that the observed / 
corresponds exactly to fj within instrumental error (see 
Eig. 3). Most importantly, because we make appropriate 
use of the collecting lenses, the behavior of the bolometer 
signal is commensurate with the emitted integral power 
although we neglect the light loss owing to the Eresnel 
reflections at the lens surfaces. 

Eigure 4 shows the emission power from the outermost 
IVC curve versus V at = 15, 20, 25, and 30 K at differ¬ 
ent offsets. The inset in Eig. 4 shows the dependence 
of the outermost IVC curve. With increasing T^, the 
IVC hysteresis and the concomitant maximum voltage 
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FIG. 3. Emission frequency measured using the interfer¬ 
ometer system (left) and the detected emission power (right) 
versus voltage. The dashed line represents the Josephson re¬ 
lation f — fj — {2elh)VlN for N = 851. 


are suppressed because of the reduction in the quasipar¬ 
ticle tunnel resistance. The voltage condition for the TM 
(m,p) mode can be obtained using Eqs. (1) and (2) with 
fj = fmp 


^ mp 
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The six vertical lines represent the calculated V^p for 
(0,2), (1,2), (2,0), (2,1), (1,3), and (2,2), respectively. 
Note that in the calculation of V^p^ we use only the realis¬ 
tic sample dimensions (re, £) and e = 17.^1^. We observe 
that the emission power has local maximum at each V^ p 
and its peak height varies depending on The slight 

difference between calculated V^ p and observed peaks of 
emission power arises from the typical 5-10% difference 
in the top and bottom lengths of the actual stack due to 
the trapezoidal cross section. Particularly, the TM (2, 0) 
and TM (2,1) modes have relatively high peaks in these 
modes. This observed behavior is consistent with the 
simulation results, which is discussed later. Note that 
the TM (1,0) mode cannot be identified in Fig. 4, since 
the resonance voltage of V^q = 0.631 V falls short of the 
applicable voltage range. 

Each emission peak in Fig. 4 has a nearly symmet¬ 
ric width of AV^p that may reflect the quality of the 
cavity resonance. None of the obtained AV^p values 
is smaller than the calculation error of V^p that stems 
from difference in the top and bottom length of the ac¬ 
tual stack. An effective Q-value can be estimated from 
Qmp = example, = 27 at = 30 K. 

The observed low-Q cavity resonance suggests that the 
broadening effect of the emission is independent of the 
angle of the sidewall inclination in the mesa but is caused 
by the impedance mismatching. Moreover, a mixing ex¬ 
periment showed a 23 MHz spectral linewidth at the min- 
imunJi^^ With careful consideration of the minimum res¬ 
olution limit of 0.75 GHz of the FT-IR spectrometeJl^, 
a Q-value in the order of 10^-10^ is estimated, which is 
much larger than Q^p. This high-Q resonance is inter¬ 
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FIG. 4. Observed emission power from the outermost IVG 
curve at different offsets versus voltage at various T^. Dashed 
lines represent the calculated using Eq. (3) for the TM 
{m,p) modes. The inset shows the Th dependence of the out¬ 
ermost IVG curve. 


preted as a unique feature of the synchronized IJJ stack; 
however, further synchronization studies are needed to 
provide supporting information. 

The Joule heating over the entire mesa comes off the 
backbending IVC curve in the resistive state. In the high- 
bias regime, the temperature distribution in the emitting 
stack is inhomogeneous due to the enormous local heat- 
in^m^ This would lead to a shorten length of the inter¬ 
nal cavit}im. Nevertheless, the overall peak structure in 
Fig. 4 is insensitive to T^, which implies that the effective 
cavity condition including re, and e is nearly indepen¬ 
dent of Ti) in the measurement range. 


IV. SIMULATION 

We use the commercial simulator SonnelP^ to analyze 
the experimental result in terms of the cavity condi¬ 
tion. By using the method of moments with Sonnet, we 
can calculate the EM field from a numerical integral of 
Green’s function at the discretized surface of the patch. 
The current density at the patch surface can also be 
calculated by substituting the EM field into Maxwell’s 
equations. Although the simulation results hereinafter 
presented do not allow us to take into account the syn¬ 
chronization among the numerous IJJs, which is actually 
very important for the physics of coherent emission, we 
find that they are consistent with the experimental re¬ 
sults and thus informative to readers. 

The simulated objects consist of the infinitely-thin su¬ 
perconducting patch atop the isotropic dielectric layer 
with thickness t. We use the dielectric constant of 
e = 17.6 from Refill. The lateral sizes of the cell 
and the patch in the analytical space are 5x5 and 
100x140 /im^, respectively. The patch is surrounded by 
electrical grounding 5 pm apart from the edges, which 
is consistent with the actual situation (see the image in 
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FIG. 5. (a) Simulation of the reflection coefficient of a 100xl40-)am^ mesa, obtained from the scattering analysis, (b)-(g) 

Snapshots of the amplitude distribution of the in-plane surface current density Jxy at the resonance frequencies. Dashed lines 
in each panel represent positions of the line nodes of the EM standing wave. 


the inset of Fig. 2(b)). The EM feeding point with an 
area of S' = 5 x 5 fivo? is set near the center of the patch. 
Note that the feed position has insignificant effect on the 
calculation result of the resonance frequency but affects 
the quality of the resonance according to the impedance 
matching. 

Figure 5(a) shows the reflection coefficient at the feed 
point in the range of 0.4-1.0 THz. Sharp drops are 
clearly observed at 0.516 THz, 0.638 THz, 0.716 THz, 
0.766 THz, 0.864 THz, and 0.894 THz. These drops are 
attributed to the EM emissions from the mesa to the 
free space. The simulated resonance frequencies are con¬ 
sistent with calculated from Eq. (2). In Fig. 5(a), the 
reflection coefficient is normalized by the characteristics 
impedance Z at the feed point, which is determined in 
general by the ratio of the voltage to the current. We 
assume Z = vj where v — VrjN is the voltage per 
IJJ and Jc = Iclw£ is the amplitude of the AC Joseph- 
son current. Using realistic parameters Vr = 1.0 [V], 
N = 850, and Ic = 18 [mA], we obtain Z = 36.7 U. In 
this particular case, we attribute the most intensive emis¬ 
sion at the TM (2,0) mode to the impedance matching. 

Snapshots of the amplitude distribution of the surface 
current density J^y at the resonance frequencies are dis¬ 
played in Figs. 5(b)-(g). The surface current is enhanced 
at the edges of the mesa owing to the edge singularity. 
The simulated J^y is associated with the standing wave 
inside the mesa: Jxy = 0 and dxyJxy = 0 correspond 
to the anti-node and the node, respectively. The dashed 
lines in each panel represent line nodes of the standing 
wave. We can identify the excited TM (m,p) mode from 
the numbers of the line nodes as follows: (b) (0,2), (c) 
(1,2), (d) (2,0), (e) (2,1), (f) (1,3), and (g) (2,2). 


V. CONCLUSION 

The nearly square Bi-2212 mesa emits coherent and 
broadly tunable CW terahertz waves depending on the 


internal EM cavity resonance. We used a collection 
optical system to directly identify the excited cavity 
modes. The /(U) that was simultaneously measured us¬ 
ing a wedge-type interferometer verified the validity of 
the Josephson relation. We analyzed the experimental 
data using the field simulator, and found that EM stand¬ 
ing waves could form along the width and length of the 
mesa. The most important new aspect of using a nearly 
square stack is that we can control the excited cavity 
mode and the concomitant polarization using an identical 
source. Thus, compact and solid-state terahertz sources 
using the stacks of IJJs are thought suitable for coherent 
communication applications. 
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Appendix A: Collimation effect of the silicon lens 


In order to collect the emitted integral terahertz waves 
from the sample, we used a collection optical system that 
consisted of a hemispherical silicon lens with a diame¬ 
ter of 5 mm, plastic lenses (PAX Co., Terahertz-Super 
Lens “Tsurupica”), and off-axis parabolic mirrors. Fig¬ 
ures 6(a) and 6(b) show the far-held emission patterns 
that were measured separately using another system with 
and without the silicon lens, respectively. In this case, 
the silicon lens was attached atop a conventional rectan¬ 
gular mesa with 60x300 /im^ that emits at the TM (1,0) 
mode. 

The difference between the two cases clearly proves 
the collimation effect of the silicon lens that makes the 
pattern more unidirectional. Based on simulations, the 
far-held pat tern was found to be ahected by the excited 
cavity modd^^^. The collecting angle of the plastic lens 
is more than 30°; hence, the integral emission power from 
the sample can be monitored independently of the excited 
cavity mode. 
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FIG. 6. Far-held emission pattern at the TM (1,0) mode 
measured (a) with and (b) without the collimation silicon 
lens. The detection angle varies in the E-held plane parallel 
to the mesa width. 


Appendix B: Wedge-type interferometer system 

In the measurement of the emission frequency /, we use 
a handmade wedge-type interferometer. A wedge con¬ 
sisting of two quartz plates and a stainless-steel spacer 
was hxed on the ID scanner, as shown in Fig. 1. The 
raster scan was performed in the horizontal {X) direc¬ 
tion, where the interspace distance of the wedge d in¬ 
creases linearly with X. The emitted monochromatic 
wave allows us to observe the interference pattern of the 
bolometer output. As an example, we plot the bolome¬ 
ter output versus X at 1.11 V in Fig. 7. The solid line 
represents a sinusoidal wave function at the best fitted 
frequency from the least squares analysis. Using this in¬ 
terferometer system, we can obtain / with an uncertainty 
of 2% much faster than using a commercial FT-IR spec- 
trometei^. 



X (mm) 


FIG. 7. Interference pattern obtained using the wedge-type 
interferometer. The solid line represents a sinusoidal wave 
function at the best fitted frequency of 0.633 THz. 
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